Periodic-cylinder vesicle with minimal energy 
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Abstract 

We give some details about the periodic cylindrical solution found by Zhang and Ou-Yang in 
[Phys. Rev. E 53, 4206(1996)] for the general shape equation of vesicle. Three different kinds of 
periodic cylindrical surfaces and a special closed cylindrical surface are obtained. Using the elliptic 
functions contained in mathematic, we find that this periodic shape has the minimal total energy 
for one period when the period-amplitude ratio j3 ~ 1.477, and point out that it is a discontinuous 
deformation between plane and this periodic shape. Our results also are suitable for DNA and 
multi- walled carbon nanotubes (MWNTs). 
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I. INTRODUCTION 



Lipid bilayer membranes have abundant shapes in aqueous environment. The equilibrium 
shape of a vesicle is determined by the minimization of the Helfrich free energy: M F H = 
J f[^(2H + Cq) 2 + kA]dA. Where H, A and Cq are the mean curvature, Gaussian curvature 
and spontaneous curvature, respectively, dA is the surface area element, and k c and k are 
constants. For close shapes, we need the surface constraint and volume constraint: 

Fs = A § dA and Fy = AP § dV, respectively. Here dV is the volume element for the 
vesicle, A and AP are constants. The total energy is F = Fjj + F$ + Fy. By studying 
the first variation of the total energy S^F = 0, Ou-yang and Helfrich obtained the general 
shape equation of closed membranes^ 2 ' 3 ' 

k c (2H + C )(2H 2 -2A- C H) 

—2XH + 2k c V 2 H + AP = 0, (1) 

where the V 2 is the Laplace-Beltrami operator. Eq.([T]) is a high order nonlinear partial 
differential equation. It is a hard task to solve it analytically. Only the general solution for 
cylindrical vesicles has been discussed^ 4 ' 5 ' and other special solutions including the Clifford 
torus,® the discounts,^ and the beyond-Delaunay surface® have been given. Since the 
numeric methods, such as finite element method using the Surface Evolver,^ were used to 
search solutions for this equation, many interesting shapes were obtained in Refs. [10-18]. 
Compared with those complicate shapes obtained by numeric methods, present known ana- 
lytical solutions are only the tip of iceberg in the abundant solutions of Eq.flT]). In addition, 
a generalized model with an arbitrary energy density function T = T[H^ A) was considered 
and the corresponding equilibrium equation was obtained. 1191 Moreover, in Refs. [20-22] , Tu 
and Ou-Yang obtained the basic equations for open vesicle with free edges using the exterior 
differential method. These basic equations are also complicated and we can suppose that 
they also have abundant solutions. 

II. PERIODIC-CYLINDER VESICLE WITH MINIMAL ENERGY 

In cylindrical case, Eq.([T]) is reduced to an elliptical equation which was fully discussed 
in Ref.[5]. Let a cylindrical surface along the y axis and tan^rr) = here ip(x) is the 



2 



angle between the x axis and the tangent to the curve at point x, the mean curvature and 
Gaussian curvature of this surface are 

1 , dib 1 d(sinib) , , N 

Besides the planar solution dip/dx = 0, the shape equation is reduced to 

d 2 9 A , ,dg 2 sec 2 V> 

#^- 4tan ^#- (1 - 2tan ^=^p (3) 



where (7 = (4^) 2 and x defined by^A- = -A + Here we let x = ./ , i'- , Furthermore, 



elastic theory is also used to analysis the shape of DNA chain and MWNTs. In planar case, 
total energy of DNA' 23 ! and MWNT' 24 ' can be written as F = % f(K - C ) 2 ds + A / ds (K 
is the curvature of the central line of DNA and MWNT and ds is the element of arc length), 
which induces the same shape equation as Eq.(j3J). Zhang and Ou-Yang found a solution for 
Eq.Q with the form' 4 ' 

here a and C\ are constants and we choose C\ = because it is corresponding to a remove 
along x axis. Choosing the reduced value X = x/xq, expression (j4j) can be written as 

sin-0 = —X 2 — a. (5) 
4a 

The cross section of this cylinder can be obtained by 

Z(X)-Z(X l )=f tan^dX, (6) 

JXi 

where Z = z/xq. Choosing different a, we find that this solution can give us four kinds of 
different shapes as been shown in Fig{TJ Here we only discuss the case a > 0, because we 
will obtain the similar results when a < 0. Valid shapes are only in the range < a < 0.462. 
When a > 0.462, the shape will be self-intersected. In FigJ21 we give a new kind of periodic 
shape using FigJT](c). Specifically, when a ~ 0.652, the points A and B in FigOJa) will 
be superposed and we get a closed cylinder with the width-high ratio D/(2X m ) ~ 0.328 in 
FigETJb). Simultaneously we need 

-X m < X < X m for < a < 1, (7) 
-X m < X < X n or X n < X < X m for a > 1, (8) 



FIG. 1: Four kinds of different shapes obtained by choosing different a in expression ([5]). (a) A 
period of a kind of periodic cylinder when < a < 0.462; (b) A special closed cylinder when 
a = 0.652, the width-high ratio D/(2X m ) ~ 0.328; (c) A period for a new kind of self-intersected 
periodic cylinder when 0.652 < a < 1; (d) A period for a kind of self-intersected periodic cylinder 
when a > 1. 



FIG. 2: A kind of periodic shape (three periods) by choosing a = 0.9 in expression ([5]). This shape 
and the shape in Fig. 2 of Ref.[4] are close to the experimental shapes in Fig.l of Ref.[25]. 

where X m = 2^Ja(a + 1) and X n = 2^/ 'a(a — 1). 

Now, we turn to discuss the total energy of this solution. For simplicity, we only consider 
the case < a < 1. Then, the total energy (actually the energy density along y axis) for a 
period is 



X 




z 




(9) 
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FIG. 3: The curve of T(a). It has the minimum T m i n ~ 2.235 at a ~ 0.256. When a — > and 
a — > 1, we have T — ► oo. 



here we define 



9 = i arcsinh \ I — - — , k 



1 — a ' 



a + 1 



the function T(a) is 
T(a) 



= {(a + l) x E[0,k] - (2a + 1) x F[0,k]}, 



(10) 



(11) 



A/a(a + 1) 

where F[x,y] and E[x,y] are the elliptic integral of the first kind and the seconde kind, 
respectively. We find that T(a) has the minimum T min ~ 2.235 at a ~ 0.256 (see Figj3]). 
Thus, the minimal total energy for one period is 



8.94 A /A; C 2 C , 2 + 2Xk c . 



(12) 



Sense it is difficult to discuss the stability of this infinite periodic solution, we only can think 
that the stable shapes satisfy a ~ 0.256. 
The reduced period T is 

8ia 



T(a) = 2 / tan^ciX 

-x m A/«(a + l) 



[a + 1) x E[0,k] -F[6,k) 



(13) 



We show T(a) in FigJH and find that T(a) has the maximum T max ~ 1.672 at a ~ 0.256. 
Here we should note that the real period which can be measured in experiment is T' = xqT 
and we have T' max ~ 1.672xo- Further, it is interesting to know why T\{a) and T(a) have 
the minimum and maximum when a ~ 0.256, respectively. However, because the elliptic 
function is everywhere continuous but everywhere non- derivable, it is difficult to discuss this 
coincidence. 
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FIG. 4: The curve of reduced period T{a). It has the maximum (for not self-intersected shapes) 
T max ~ 1.672 at a ~ 0.256. When a -> 0, we have T -> 0. 

An important characteristic parameter for this periodic shape is the period-amplitude 
ratio, it is 

/3(c) = -ji^— (« + l)xE[^k]-F[fl,k] . (14) 

We find that f3{a) has the maximum (for not self-intersected shapes) f3 max ~ 2.4 when a — ► 
(see Figj5j), which is in accord with the value in Ref. [4]. Specially, when a = 0.256, we get 
(3 = 1.477. Through measuring the shapes in Ref. [26], we obtain /3 ~ 2 for Fig. 9, /? ~ 3.3 for 
Fig. 11(a), and (3 ~ 2.7 for Fig.ll(b) and (c). It indicates that not all these shapes can be 
explained by this solution. Moreover, in Ref. [26] the authors found that planar surface can 
gradually change into periodic cylinder. But in this solution, it seems that the planar surface 
cannot change into periodic cylinder continuously, because planar surface means (3 — * oo 
which is out of the range (3 < 2.4. Considering planar surface is always stable scampered 
with any periodic surface, ' 27 ' we think there is an energy barrier between this two kinds of 
surfaces. Supposing that a plane with Cq and A changes into a periodic cylinder with Cq 
and A, then the energy change between this deformation can be obtained by the following 
way. 

For this periodic cylinder, the length for one period of the cross section line is 

Lx(a) = 2x [ -^—dX = -i 8x J^—x F[0,kl. (15) 
i-x m cos^ \ a + 1 

Consider the energy density for planar surface solution is T v = yC*o + K the energy change 
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FIG. 5: The curve of period-amplitude ratio (5{a). It has the maximum (for not self-intersected 
shapes) (3 max ~ 2.4 when a -> 0, and 0(0.256) ~ 1.477. 

for a period is 

AT V = Tiia) - FpL^a) 

= ^(a)-(^ + A)Lx(a). (16) 

Let Co = Co and A = A, there is 

AJp = 4^ c 2 C 2 + 2AA; c ft(a), 

with 

fi(a) =r(a) + i — x F[0,kl. (17) 
V a + 1 

We show Q(a) in Figj6] and find that Q(a) has a minimum f2 m j„ ~ 1.325 at a ~ 0.652, 
which is corresponding to the closed shape in Fig. 1(b). Simply let the barrier be equal to 
the energy difference, then the minima barrier is between the planar and the self-intersected 
8-shape in Fig{T](b). Clearly, in phase space, if we want that a planer drop to a periodic 
cylinder, we need lS.T v < in Eq. flTBT) . It cannot be satisfied if Co = Co and A = A, because 
Q(a) > 0. Thus, in this downward transition process, A (or Co) will change. 

When a > 1 we find that .Fi(a) is monotone decreasing following the increase of a and 
T\ = 4nk c Co when a — > oo. Another way, when AP = 0, a circular cylinder solution 
with radii R needs C R = l. [28] Its total energy is T = §^{1/R+ C )ds = 47rA; c C which 
is in accordance with the above result (note that a — > oo is corresponding to a circular 
cylinder! 4 !). 



7 



Q(a) 
5 

4 . 5 

4 

3.5 

3 

2.5 

2 

1 . 5 

0.2 0.4 0.6 0.8 1 

FIG. 6: The curve of 0(a). It has the minimum £l m in — 1-325 when a ~ 0.652. When a — > 
we have — > oo and when a -» 1 we get O ~ \/2- In addition, we have 0(0.256) ~ 1.526, 
0(0.462) ~ 1.356. 

Remarkably, although the shape in FigOJb) is self-intersected for 2D surface, it is useful 
to describe the DNA and MWNT shapes, because that any perturbations on the binormal 
direction will probably lead it to be not self- intersected. A ring solution: x = Rsinip for 
Eq.(j3J) induces 

XR 2 + h c (C 2 R 2 - 1) = 0. (18) 

If a DNA ring changes into the shape in FigETJb) and Co and A change to Co and A, 
respectively, we have 

2^ = ^(0.652). (19) 
Eqs. flIB"]) and ( fl~9"l) yield the relationship 

2.846 x (C 2 + 2A/A; C ) = C 2 + 2A/A; C . (20) 
If Co = Co = 0, expression (120"]) indicates that the tension A will increase in this deformation. 

III. DISCUSSION AND CONCLUSION 

In conclusion, we have studied the particulars about the periodic cylinder in Ref. [4] for 
closed shape equation of vesicle. Several characteristics for this periodic shape are obtained, 
which give us some interesting details about this solution. For instance, when the period- 
amplitude ratio (3 ~ 1.477, this periodic shape has minimal total energy for one period. We 
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hope the shape with period- amplitude ratio (3 ~ 1.477 can be found in experiment. As for 
the other periodic surface: beyond-Delaunay surface J 8 ' because there isn't analytic integral 
function of its total energy, we don't know whether there is similar result. 

Finally, we point out an interesting revelation. We know that if the energy density of an 
one-dimensional structure, such as DNAt 23 ] and MWNTs,' 24 ' can be written as 

F= l -k c (K-C Q )\ (21) 

the corresponding shape equation in planar case is 

2K + K 3 - (Cl + 2\/k c )K = 0, (22) 

where K = dK/ds. This equation is equal to Eq.Q (except the straight line solution 
K — 0). The shapes of DNA and MWNT in planar case satisfy this equation. Thus, the 
solutions of Eq.([3]) have broad meaning. Ou-Yang et al. proved that the energy density of 
MWNTs can be written as T = a if 2 J 24 ' but this model seems not suitable for single- walled 
carbon nanotubes (SWNTs), because the constant a will be zero for SWNTs. However, we 
note that the shape in Fig. 2 of Ref. [4] and the new shape in Figj2]of this paper seems to be in 
good agreement with the experimental SWNT shapes in Fig.l of Ref. [25]. Thus, we suppose 
that the energy density of SWNTs satisfies expression (1211) under suitable approximation. 
We will try to prove this supposition in our future work. 

This work is partly supported by National Science Foundation for Young Scientists of 
China, Grant No. 10804129. 
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